On Electron-Positron Pair Production by a Spatially 

Nonuniform Electric Field 

A. Chervyakov * and H. Kleinert ' 
Institut fiir Theoretische Physik, 
Freie Universitat Berlin, 
Arnimallee 14, D- 14195 Berlin 



Abstract 

A detailed analysis of electron-positron pair creation induced by a spatially nonuniform 
and static electric field from vacuum is presented. A typical example is provided by the 
Sauter potential. For this potential, we derive the analytic expressions for vacuum decay 
and pair production rate accounted for the entire range of spatial variations. In the limit 
of a sharp step, we recover the divergent result due to the singular electric field at the 
origin. The limit of a constant field reproduces the classical result of Euler, Heisenberg 
and Schwinger, if the latter is properly averaged over the width of a spatial variation. The 
pair production by the Sauter potential is described for different physical regimes from 
weak to strong fields. In all these regimes, the locally constant-field rate is shown to be 
the upper limit. 

1 Introduction 

Since Klein's formulation [I] of his gedanken experiment in 1929 and its interpretation by 
Sauter [2], Heisenberg and Euler [3], and Hund p_] it became clear that the vacuum of quantum 
electrodynamics (QED) is unstable against the electron-positron pair production in the present 
of external electromagnetic field. The first discussions were all done at the level of first-quantized 
field theory. The results were rederived within second-quantized field theory in the one-loop 
approximation by Schwinger [5] for a constant electromagnetic field, and Nikishov [TJ [8] and 
Narozhny and Nikishov [9j [[_)] for more general field configurations. The two-loop radiative 
corrections were calculated by Ritus [HI [12] and Lebedev and Ritus [13]. The spin-statistics 
connection in QED with unstable vacuum was discussed by Feynman [13]. A general discussion, 
detailed calculations and further references on this subject can be found in Refs. [151 HHJ El 
HH [TH [201 EH [221 123]. For recent developments and the modern state of the problem, see 
Refs. [2UE5]. 
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In a constant electric field E, the probability for vacuum decay via pair creation in unit 
volume per unit time (i.e., the vacuum decay rate) can be written as 



where E c = m 2 c 3 /\e\h ~ 1.3 x 10 18 V/m is the critical field, whose work to move an electron over 
the Compton wavelengthA e = h/mc being equal to its rest energy mc 2 andAg/c ~ 7.3xlCT 59 m 3 s 
is the Compton space-time volume. The rate ([I]) results from the imaginary part of the one-loop 
Euler-Heisenberg Lagrangian [3] . The formula (pQ) was anticipated by Sauter [2J and derived by 
Euler and Heisenberg j3] and, more elegantly, by Schwinger [5J. It is presently often referred to 
as Schwinger formula, and the pair creation process as Schwinger mechanism of pair production. 
The process can be understood within the Dirac picture of vacuum as a quantum tunneling 
of electron through an electrostatic potential barrier created by the electric field. Its rate is, 
however, exponentially suppressed for presently accessible field strengths E achieved roughly 
4 x 10 14 V/m by optical lasers [26J. This leads to the present impossibility of observing the pair 
creation for which the extraordinary strong electric field strengths E of the order or above the 
critical value £J C ~ 1.3 x 10 18 V/m are required. 

While creating such strong macroscopic fields in the laboratory is impossible today, it will 
hopefully be achieved in the near future after further advances in the laser technology. Indeed, 
the energy extraction and the optical focusing can be improved considerably in the X-ray free 
electron lasers to produce the electric field of the order of E ~ 10~ 2 E c capable of yielding a 
sizeable effect of pair production [271 12H] • Under these circumstances, we could get a significant 
progress in understanding the properties of QED beyond the scope of perturbation theory in 
the strong-field regime. Of course, such strong macroscopic fields in the laboratory would never 
be as uniform as those originally considered in the Schwinger mechanism of pair creation. Their 
spatial and temporal modulations must be taken into account in optical and, especially, in X-ray 
free electron lasers (for review, see Ref. [26] and for recent estimates of pair production created 
by laser fields, Refs. [291 (30J [311 [321 [331 El])- Apart from these, enormous electromagnetic and 
gravitational fields are expected in the powerful Gamma Ray Bursts, strong Coulomb fields in 
the relativistic heavy-ion collisions and nonabelian gauge fields in QCD, for which the Schwinger 
mechanism of pair creation is also relevant [23]. 

Going beyond the constant fields in space and time is therefore an important goal in studying 
pair creation. First estimates of the impact of spatial and temporal inhomogeneities on pair 
production can be found in Refs. [HI [7J for spatially varying fields, in Refs. [H QUI ESI ESI S3 
[38[ [39] for oscillating electric fields, and in Ref. [30] for an electric field enclosed by conducting 
plates. Recent studies of pair production by inhomogeneous electromagnetic fields including 
developments of many useful approximate schemes such as WKB method, world-line formalism, 
instanton and Monte Carlo techniques can be found in Refs. [3111321,1331133113511151137113811391 
1501 I5T1 [521 [531 153] . The back-reaction of produced pairs on the external field has been also 
taken into account in Refs. [55J, |56] , and more recently in Ref. [57] . 

For an arbitrary shape of background field the study of pair creation becomes rather in- 
volved, so that only special field configurations have been considered. In fact, one has to 
find the scattering and pair production probability densities in terms of solutions of the single- 
particle Klein-Gordon or Dirac equation in the presence of external field yielding the Bogoliubov 
transformation coefficients within the S'-matrix formalism. In addition, the energy-momentum 
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integral over the level-crossing (Klein) region involving logarithm of the reflection coefficient 
must be performed. In this way, one can connect in a most transparent and efficient way the 
one-particle Dirac theory, in which the Klein paradox appeared with the second-quantized field 
theory, in which it was resolved satisfactorily 

In this paper, we examine the problem of pair creation of Dirac particles in a spatially 
nonuniform and static electric field E(z) = (0,0, E(z)) directed in the z-axis and varied only 
along this axis, where E(z) is the Sauter field of the form 



with the peak E = vk/\e\ at z = 0. The parameters v and 1/k are associated with the height 
and the width of the electrostatic potential V(z) = e<fi(z) created by the field E(z) = —(f)'(z), 
respectively. Explicitly, it reads 



with v, k > 0. The limit k — > oo leads to a sharp step, whereas the limit k — > with vk fixed 
reproduces the linear potential caused by a constant field E . 

With no magnetic field, we use the gauge A^(t, xj_, z) = (4>(z), 0, 0, 0) in which the Dirac 
equation allows the full separation of variables. As a result, the problem of pair production 
is reduced to an exactly solvable quantum-mechanical problem [HI CO E] of the scattering of a 
Dirac particle on the static potential barrier ([3]). The problem is recapitulated in Section 2. 
This provides us with exact expressions for reflection and transmission coefficients in terms of 
scattering data for both the sub- and the supercritical values of the potential fl3]). For super- 
critical barriers with the level crossing between positive and negative energy levels, the Klein 
paradox is connected with the vacuum instability via the electron-positron pair production. In 
Section 3, we discuss how the reflection and transmission amplitudes are related to the pair 
production rates, first for a certain particular state, and then by taking all these states into 
account. In particular, the total vacuum decay rate due to pair creation is expressed as an 
energy and transverse momentum integral over the level-crossing (Klein) region involving the 
logarithm of reflection coefficient. 

The quantum-mechanical scattering amplitudes on the potential barrier fl3]) have been known 
since Suater's work [2], and their relation to the pair production probabilities since Nikishov's 
work [SI El [H]- In order to find the pair production rates, we must however complete that work 
by evaluating the energy-momentum integral for a given local probabilities. The method of 
calculating such an integral has been developed recently in Ref. [53]. In Section 3, we apply 
this method to perform the rotationally invariant integral over transverse momenta in three 
dimensions exactly. The remaining integral over the energy, after an appropriate change of 
variable, provides us with simple spectral formulas for the vacuum decay rate and also for the 
pair production rate. The analytic expressions for these rates are derived in Section 4 for the 
entire range of the parameters v and k of the supercritical potential ([3]). This allows us to 
access different physical regimes from weak to strong fields. 

In the weak-field regime corresponding to presently available field strengths, the supercritical 
potential fl3]) can only create the pairs near the constant limit. The pair production rate is very 
small and is described by the semiclassical expressions of Refs. [ISJ H5J H7J HSJ HHJ ED] which are 
recovered and slightly corrected in this paper. The constant-field limit of the vacuum decay rate 
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is found to be different from the Schwinger formula ([I]). It agrees, however, with formula ([T]), 
if the latter is properly averaged over the width of a spatial variation. As a such, this limit 
represents the locally constant-field rate. The explanation as well as the detailed comparison 
of the two constant-field rates can be found in Appendix A. 

With increasing field strength from weak- to strong-field regime, the supercritical poten- 
tial ((3]) creates the pairs much more intensively over entire field width extended from sharp 
to the constant limit. The corresponding expression for the vacuum decay rate becomes very 
large but interpolates analytically between these limits always below the locally constant-field 
rate. The latter is therefore the upper limit for the spatially nonuniform electric field (j2]). 



2 Barrier scattering and pair production 

Consider a relativistic spin- 1/2 particle of charge e (e = — |e| for electron), mass m and energy e 
moving along the z-axis in the field (T5]) of the potential barrier ([3]). The energy and momentum 
of a particle in this potential read (c = 1): 



Po(z) = £ - V(z), p 3 {z) = \bl{z) - (p 2 ± + m 2 ), p ± = Pl +p 2 



(4) 



In the transverse direction the particle propagates freely as a plane wave exp[2(p_i_xj_ — et)/h}. 
Its wave function can therefore be separated as a product of this plane wave with a four- 
component spinor field ip(z) satisfying the one-dimensional Dirac equation 

[ 7 °Po(z) + 7^P± + ^d z - m] il>{z) = , (5) 

where 7 M are the Dirac matrices. The solution of Eq. (jSJ) is represented in the following form 

%j)(z) = [l°Po(z) + 7 ± p_ L + ih^d z + m] x(z) , (6) 

where the four-component spinor field x( z ) satisfies the second-order Dirac equation 



This has two sets of solutions, each associated with two spin states a = 1, 2 of a particle 

ip + i(z)u a , a = 1,2 



(7) 



(p-i(z)v a , a = 1,2 



(8) 



where the constant spinors u a and v a are eigenvectors of the matrix (7°7 3 ) with eigenvalues +1 
and —1, respectively. In the Dirac representation of 7-matrices, these read explicitly, 
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(9) 



with the normalization u\u a i = 8 aa >, v\vj = 8 aa > and u\v a i = 0, v\u a i = 0. 
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The functions <p±x(z) in Eq. (jSJ) obey the Schrodinger-like equations 



<p"( z ) + ^ IpK z ) - ™heE(z)] ip s {z) = , s = ±1 . (10) 

For each s, there are two independent solutions which behave asymptotically as plane waves 
outside the potential barrier ([3]). Let us introduce the initial and final values of the particle 
energy far to the left and to the right of the potential 

p ( o f) = Po(z)\ z ^ 00 = e±v, (11) 

and the corresponding momenta 

4 T) = P3(*)U TOO = VPo T)2 -(pi + ™ 2 ), (12) 

where we take the positive square roots. We assume also that momenta pg in Eq. (1121) are 
real thus restricting the asymptotic energies to \p^ \ > \fp\ + m 2 . 

For both signs s = ±1, the independent solutions possess the same asymptotic behavior at 
z — > — oo as well as at z — > +oo due to vanishing the electric field (J2j) . Let ±ip s {z) be the 
two independent solutions that describe the ingoing (+) and the outgoing (— ) plane waves at 
z — > — oo, 

^WLh»«^. * = ±1 > ( 13 ) 

whereas (p a (z) the two independent solutions that describe the ingoing (+) and the outgoing 
(— ) plane waves at z — )■ +oo, 

^WL+^e^V*, s = ±1 . (14) 

Depending on two projections r = ± of momenta p^ onto the z-axis, we obtain two complete 
sets of solutions { r (p s } and { r f s } for equations ( fTUl) satisfying the asymptotic conditions ( fT3l) 
and (TH|) . respectively. By virtue of linearity of equations (flOj) . each function from the one set 
{±<p s } can be expressed as a linear combination of the both functions from the another set 
{^tfis} and vice versa. The coefficients of these combinations are then related to the reflection 
and transmission amplitudes of a particle in a scattering problem. Explicitly, these can be 
found for potentials for which the equations (1101) are solvable exactly. 

For each s, the functions {±tp s } as well as the functions {^tps} are linearly independent. For 
each r, the functions { r (p±i} as well as the functions { r tp±i} are related to each other. The 
relations for r (p + i(z) and r {p_i(z) follow from Eqs. (fTOl) and (Tl4|) . Explicitly, these read 

[m z +p (z)\ Vi(z) = [pI +) TP { 3 +) ] Vi(*), (15) 

w-^ivw = -U +) ±ri +) l ww. (is) 



Similar relations can be found for r (p + i(z) and r </?_i(;z) with the help of Eqs. ( TlOl and ([131) . 

According to Eq. (IE]), each complete set of solutions { r (p s } and { T Vs} provides eight solutions 
for the second-order equation (J7j) with only four solutions being linearly independent. Using 
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for instance { r <p s }, we obtain eight solutions explicitly as r ip + i(z) u„ and r <£>__i(z) v a with r = ± 
and a = 1,2, where the functions r ip + i(z) and r {p_i(z) obey Eqs. f fT5|) and ffT6l) . Thus, we can 
use either r <y9 +1 (z) u^, or r <^_ 1 (z) v a in order to obtain four independent solutions for the Dirac 
equation ^ with the help of Eq. ([H]). Substituting for instance ± (p^.i(z)v a into Eq. ^ and 
applying Eqs. ffT5l) and (|T6l) yields one complete set of solutions as 



± tfj a (z) = [j°p (z) + 7 X P± + rtl 3 d z + m] ± ifi-i(z) v c 

= ± n{ ± ^ 1 (z)^+ [pi +) ±4 +) ]ViW^} , <r = l,2, (17) 

where £ CT = (7 _L pj_ + 77i)i><7 and ?7 CT = (7°v (J ) for a = 1, 2 and ^iV are normalization constants to 
be determined later. Similarly, the another complete set of solutions for the Dirac equation (j5j) 
reads 

±ip a (z) = [7% (2) + 7 ± p_ L + ihfd z + m] ±¥>-iCz) «<r 

= ±iv{±^_ 1 (z)^+ [p^ip^JiV+i^)^} , a = 1,2, (18) 

with the normalization constants ±N. In Eqs. (ITTj) and ({TBI) , the spinor fields with a ^ a' are 
orthogonal ^ a {z)^) a i{z) = 0. 

In order to solve equations flTUJ) with the potential (j3J), we set 

/4 _) = 2/ifc/z , j4 +) = 2Hkv , (19) 

where 



" 2 - 2 = A sp < 20 > 



with 



A,^>0. (2!) 



In these notations, equations ( JTU|) become 

^±i (*) + — Si- [" 2 e2fcz + ^ e ~ 2fc * + [y 2 + - ^ 2 ) ± «A] (p±i{z) = o . (22) 

cosh kz 

Thus, the solution <f+i(z) = <£>(z, —A) can be obtained from the solution (f-i(z) = y?(z, A) via 
the interchanging A — > —A with A > 0. The function (p(z, A) satisfies 

<p"(z, A) + k l \v 2 e 2kz + /i 2 e" 2 ^ + (z/ 2 + fi 2 - A 2 ) - zAl <p(z, A) = . (23) 
cosh 

Equation (|23|) is solvable exactly in terms of the hypergeometric function [7] . For this equa- 
tion, we are looking for two independent solutions + f(z, A) and ~^>(z, A) with fixed asymptotic 
behaviors at z — » +oo: 

^(z, A) « exp (±2ifcz/z) , 2; -)■ +00 , (24) 
where the normalization can be neglected. Let us introduce the new dimensionless variable 

C = -exp(-2A;z), (25) 
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running from — oo to 0. To find + (f(z, A), we use the substitution 

MC, A) = (-C)-* /(C, A) , (26) 
where the function f((,\) satisfies 

U" (C, A) + (i - 2«o/'(c, A) + [^r^r " TT^r] /(c ' A) = • (27) 

To remove the singularity at £ = 1 in Eq. fl27|) . we replace 

/(C,A) = (l-C) iA MC,A). (28) 
Eq. ( 12 7p becomes now a hypergeometric equation 

C(l - CK(C, A) + {(1 - 2iz/) - [(1 - 2*V) + 2zA] C} u/(C, A) + [(A - vf - ^ w(C, A) = 0, (29) 

where the function w((, A) tends to a constant for £ — )■ (z — > +oo). With this condition the 
solution is the hypergeometric function 

w(C, A) = F [i (A - v - y) , z (A - v + fj) , 1 - 2iv ; C ] (30) 

up to some normalization factor. 

The solution + <£>(.2, A) of Eq. f l23|) reads finally, 

+¥>(*, A) = (-C)- w (l-C) iA i 71 ^(A-^-/x)^(A-^ + /U ),l-2^;C], 

C = - exp(-2fe) . (31) 

For z —> +oo (£ — >• 0), the function + (p(z, A) satisfies the asymptotic condition in Eq. fl2l]). To 
find its the asymptotic behavior for z — > — oo (£ — > — oo), we use the Kummer transformation 
of the hypergeometric function. This yields the superposition 



~(p(z, A) = a(A) +(p(z, A) + 6(A) -(p(z, A) 



with the amplitudes 



a(A) 
6(A) 



r(i - 2?v)r(-2?» 



r[i(A-i/-^)]r[l-i(A + i/ + ^)] 

r(i - 2zV)r(2z» 

r[z(A-z/ + /i)]r[l-z(A + z/-/x)] 



Here the functions +ip(z, A) and _</ , (- 2 ) A) are two other independent solutions of Eq. 
fixed asymptotic behaviors at z — > — oo: 



±ip(z, A) ~ exp (±2ikfiz) 



z — > — oo . 



(32) 

(33) 
(34) 
with 

(35) 



Their explicit forms will not be necessary further. To get the asymptotic behavior of the 
function + ip(z, A) at z — > — oo from Eq. ( |32l) . we apply the conditions (135]) and obtain 



~<p(z, A) ~ a(A) exp (+2ikfiz) + 6(A) exp (—2ik/iz) , z — >■ — oo . 



(36) 
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In order to find the second solution ~tp(z, A), we use the invariance of Eq. (123]) under the 
transformation involving a combination of complex conjugation with the substitution A — > —A 
for A > 0. Applying this transformation to the solution in Eq. (j3~Tj) yields 

-<p(z,\) = (-0 iv (l-(y x F[i(\ + v + ri,i(\ + p-v),l + 2iis;(}, 

C = -exp(-2kz). (37) 

For z +oo, the solution ~^p(z, A) satisfies the asymptotic condition in Eq. fl24|) . For z — > — oo, 
we use the Kummer transformation of the hypergeometric function in Eq. ([3j 



-<p(z, A) = b*(-\) +( p(*> A ) + A ) > (38) 

and the asymptotic conditions fl35|) . This yields 

A) « 6*(-A) exp (+2ikfiz) + a*(-\) exp (-2ik/iz) , z -)> -oo . (39) 

The equation (I22j) for <£> + i(z) = —A) can be obtained from Eq. ( 1231) via the interchanging 
A —7- —A for A > 0. The corresponding solutions + (p(z, —A) and ~<p(z, —A) are given by Eqs. (l3Tj) 
and ( 157)) respectively, with the same interchanging. For z — > +oo, their asymptotic behaviors 
are coincident with the asymptotic behaviors of the functions ± v 9 ( z ; ^) m Eq. fl2] 



± <yj(z, —A) « exp (±2ikvz) , z — >■ +oo . (40) 

Note that both equations (1221) take the same independent of A form in this limit. Accordingly, 
the asymptotic behaviors of the second set of independent solutions +(f(z, —A) and ~(f(z, —A) 
coincide with that for the functions ±(f(z, A) in Eq. ( 155)) : 

±^?(/Z, —A) ~ exp (±2ikfiz) , 2; — )• —00 . (41) 

We use these conditions to find the asymptotic behavior of the functions ^(fi(z, — A) at z — > — 00 
from Eqs. (152")) and (1581) with the replacement A — > —A. This yields 

+ f{z, —A) o(— A) exp (+2ik[xz) + 6(— A) exp (—2ikfiz) , 2 — >■ —00 , (42) 
~ip(z, — A) « &*(A) exp (+2ikfxz) + a* (A) exp (— 2ik\iz) , z — >■ — 00 . (43) 

In Eqs. (02} , (S3D and Q5S]), flgj), the amplitudes can be found from Eqs. ([55]) and ([MD with 
the help of relations 

^+A) 6( _ A) = iMz^ 6(A)i (44) 

(fi + v — A) (/i - z/ + A) 

and their complex conjugated. 

Having obtained the functions ± (p(z, A) and ± (p(z, — A), we find from Eq. ( ITTj) four indepen- 
dent solutions of the Dirac equation j5J) explicitly as 



± 



^Pa(z) 



rtz,\)£ a + ri +) ±^ +) -A) tj. 



24 +) bS +) ±p| +) l 



(7 = 1,2. (45) 
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Their asymptotic behaviors are fixed at z — > +00, where according to Eqs. (J24j) and fHOl) we 
obtain 

HM - H*{z) = " XP(± ^ +) " M) t + U +) ± , * +00 , ^ = 1,2. (46) 

( + )| J + ) _L J + )\ 1 J 



2^K r '±ri" 

In turn, the asymptotic behaviors of four other independent solutions ±ipa{z) are fixed at 
z —00, where it follows from Eq. (fl8l) with the help of Eqs. (|35l) and (I4T1) that 



= ^(z) « ± #(*) = e ^^hB k a + (pj~) ± , ^ -00 , a = 1, 2 . (47) 

2j4 ) bS ) ±Ps ) l 

The wave functions in Eq. (|45p are normalized in such a way that the absolute value of current 
density along the z-axis evaluated with respect to these functions for the states with given 
a, e and p_i_ is equal to unity. By virtue of Eq. (j5J), this quantity is the ^-independent and 
can therefore be evaluated with the help of the asymptotic expressions fj4*6|) and ( 147|) . The 
normalization is thus \ ^(z)^ i/)^(z)\ = \±ip^ {z)~f 3 ±ip„ (z)\ = 1. 

By means of asymptotic expressions (|47|) . the asymptotic behaviors of the solutions ip a (z) 
in Eq. ( 141)]) at 2; — > —00 can be found from Eqs. (13"6"]) . (I4"2"j) and Eqs. ( 159"|) . (S3]), respectively. For 
+ -?/' CT (z), this yields 

+^(z)«J(A) + ^)+i2(A)_^), z->-oo, a =1,2, (48) 

where the amplitudes /(A) and -R(A) are determined independently of the spin state a = 1, 2 
with the help of relations (144]) as 

/(A)=f ^ + '' + ?| Y%A), fl(A)=f #i^4 ) 1/2 KA). (49) 



z/ |/i + i/ — A 

Thus, we find the solution + ip a (z) satisfying the asymptotic boundary conditions 

I{X) + ^{z) + R{X)^{z) , z^-00 

-M*) — ► { ■ (50) 

~ipv(z) , z -> +00 



The process in Eq. (150]) involves the electron wave of the energy e > impacting from z = — 00 
with the amplitude /(A) on the potential barrier ([3]) with v > at early times, which is partly 
reflected back to 2 = —00 with the amplitude -R(A) at late times. The energy of electron far 
to the left is therefore positive p$ — (e + v) > \/p\+ rn? > 0, and its momentum p% is real 
and positive in this region. 

In order to describe the plane wave solution involved in this process at z — +00, we have 
to specify the height v of the potential barrier ([3]). For weak (subcritical) potentials with 
v < e — v/pj_ + m 2 , this plane wave solution with the unit amplitude corresponds to a particle 
transmitted through the potential barrier at late times and moving towards z = +00 with the 
positive energy p^ = (e — v) > \/p\ + m 2 > and with the real and positive momentum pjj . 

However, for very strong (supercritical) potentials with the height v > e + \Jp\+ m 2 there 
are no longer particles in the asymptotic region of the large positive z. In this case, the 
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energy e of impinging electron lies in the interval —v + \fp\ + rn 2 < e < v — \/p\ + m 2 
henceforth called the Klein region. The energy far to the right becomes now negative p^ = 
(e — v) < —\Zp]_+ < 0, while the momentum p^ is real and, by our assumption, is still 
positive. For v > \fp\ + rn 2 > m, the positive energy levels far to the left overlaps the negative 
energy levels far to the right in the Klein region. The plane wave solution + ip^(z) with negative 
energy in Eq. (150]) represents the transmitted wave via a quantum tunneling. This right-moving 
wave propagates with the unit amplitude backwards in time and corresponds to an incoming 
antiparticle coming in from z = +00 and moving forwards in time towards the potential barrier. 
In the process ( 1501) . the incoming from the left at z = —00 particle annihilates the incoming 
from the right at z = +00 antiparticle. 

With wave functions (|47p . we obtain the incident and reflected electron currents in the 
asymptotic region z < for the process (1501 as follows 

j»5 = |/(A)| 2 ( + # 7+ t#) = IW ( + ^V+^) 2 = |/(A)| 2 z, (51) 

j^ cfl = |i?(A)| 2 (_^7-^) = |i?(A)| 2 (-^7 3 -^) 2 = -|i?(A)| 2 z, (52) 

where the superscript p denotes the particle state and, according to Eqs. (l4"9|) . (I55j) and (15111 . 

= sinh tt(A -,/-/*) smhn(\ + v + ri > 
1 v ;| sinh27r/x sinh27rz/ ~ ' v ; 

2 sinh7r(A - v + fi) sinh tt(A + v - p) > 
1 1 n sinh27r/i sinh27rz/ ~ 1 ' 

Here and below the upper/lower sign refers to (sub/super) critical potential barrier ([3]), respec- 
tively. Note that (A ± v — /i)>0 for (sub/super) critical potential, while (A ± v + \i) > for the 
both potentials. 

The current in the asymptotic region z > is evaluated with respect to wave functions (j46p . 
This yields 

j I'T' in = ( + ^7 + ^) = ( ^ = ±* , (55) 

where j£>Q is the transmitted towards z = +00 electron current and j">o is the incident from 
z = +00 positron current directed in the positive/negative ^-direction within the asymptotic 
region z > of the (sub/super) critical potential barrier (J3J), respectively. The superscript n 
denotes the antiparticle state. According to Eqs. (jSTl) . fl52|) and fl55|) . the conservation of current 
along the z-axis reads 

|/(A)| 2 -|i?(A)| 2 = ±l, (56) 
with |/(A)| 2 and |i?(A)| 2 are given in Eqs. pi} and fl54"j) . 
The reflection coefficient is then defined by the ratio 

^ \j"z<o a \ l-^(A)| 2 sinh 7r(A — v + n) sinh 7r(A + i> — /i) 



LSI l J ( A )l 2 sinh tt(A - z/ -/x) sinh 7r(A + z/ + /i) 



> . (57) 



It is the relative probability for elastic scattering of electron by (sub/super) critical potential 
barrier (j3j. The transmission coefficient reads 

, = b?>o /n ' m | = 1 = sinh 2tt/x sinh 2nv 

~ \jz<o\ " l J ( A )l 2 sinh tt(A - -//) sinh tt(A + i/ + //) - ' 1 ' 
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For supercritical potentials, the probability fl58|) is a relative probability for incoming electron 
to annihilate the incoming positron. The coefficients fl57|) and f )58|) were first obtained by 
Sauter [2]. 

From Eqs. f )56|) . (1571) and fl58|) we find the relation between reflection and transmission 
coefficients for (sub/super) critical potential 

r ± t = 1 . (59) 

For supercritical potentials, the reflected electron current is therefore larger than the incoming 
electron current, and even without the incoming electron, the reflected current is exactly equal 
to the transmitted current. This is known as the Klein paradox [1] for which the single-particle 
description becomes inadequate, since the time evolution operator of a single-electron wave 
function is no longer unitary as far as r > 1. Contrary, within the second quantization the 
scattering matrix operator is always unitary and the Klein paradox implies that the negative 
energy continuum of the Dirac vacuum contributes to the both currents via production of 
electron-positron pairs under the influence of supercritical potential. The probability (158]) is 
then the relative probability for spontaneous production of a single electron-positron pair. 

To show this, we go back to our discussion above Eq. ( 148]) in order to find the asymptotic 
behavior of the second independent solution ~ip a (z) in Eq. ( 145]) at z — > — oo in terms of 
asymptotic expressions (]4"T]) . From now on, we shall only consider the supercritical potentials. 
In this case, applying Eqs. ( 159"]) . (T4"3"]) to Eq. ( 14"5]) . we obtain 

-^{z)^-R*{\)+rt{z)-P{\)^{z)i z^-oo, a = 1,2, (60) 

where the amplitudes I*(A) and R*(X) are complex conjugated of the relations ( 14"9"]) . 

Thus, the solution ~ip a {z) has the following asymptotic behaviors at both infinities 

( -R*(X) + rt(z)-I*(\)^(z) , z->-oc 
— ► \ ■ (61) 

The process ( 16 ip describes directly the steady state of the electron-positron pair production 
by means of amplitudes of the previous process (150]) : within the asymptotic region z < 
of the supercritical potential barrier (j3]) there are the ingoing electron wave + ifj^(z) with the 
amplitude —R*(\) and the outgoing electron wave -tp^(z) with the amplitude —J* (A). This 
wave is accompanied by the outgoing positron wave ~ip^(z) with the unit amplitude in the 
asymptotic region z > 0. The corresponding currents far to the left are 

j£S = |i?(A)| 2 ( + # 7 +V£) = l^(A)| 2 (+^7 3 +^) z = |iWz, (62) 

j^ ut = |J(A)| 8 (_^ 7 -V'3 = |/(A)| 2 (-^7 S -V'38 = -|/(A)| a *, (63) 
and, far to the right is 

j *>7 = tVvV^) = 0W-v£) 2 = z > (64) 

with the current conservation being the same as in Eq. (156"]) for supercritical potentials. 

From Eqs. ( 16~4"]) and ( 16"3"]) . the relative probability for production of a single electron-positron 
pair is 

|j*>o Ut | = 1 = sinh27r/x sinh27rz/ > 

|j, P <o Ut | l J ( A )l 2 sinh7r(A-z/- / u) sinh7r(A + z/ + /i) - ' 1 ' 
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and therefore w = t with t being the transmission coefficient for supercritical potentials found 
before in Eq. ( 155|) . 

The average rate of produced electron-positron pairs in the process (16 ip is then derived from 
Eqs. (ED and (JS2D as 

- = |j z n >°o Ut | = 1 = sinh27r/x sinh27rz/ > 

liz<ol |-R(A)| 2 sinh7r(A-z/ + /i) sinh7r(A + z/- / u) ~ 1 J 

and, together with Eqs. (1571) and (158]) . reads 

n = t/r . (67) 

It is the absolute probability for production of a single electron-positron pair in a given state 
for which t is the relative probability and 1/r is the probability that no pairs are created in this 
state. This probability is called the vacuum-to-vacuum probability and is determined in such 
a way that the total probability for production of zero and one pairs is unity (l/r)(l + t) = 1 
which follows from Eq. (j59P for supercritical potentials. Thus creating more than one pair is 
blocked by Eq. (loT?P in the agreement with the Pauli principle for fermions. Note, finally, that 
1/r is the reflection coefficient and t/r is the transmission coefficient for the process (I6Tj) . 



3 Vacuum decay and pair production rates 

In the previous section, the particle-antiparticle pair production has been described for each 
particular state n = (p, a) which is completely specified by the total energy and transverse 
momentum p = {e, p_i_}, and also by a certain spin projection a. For these conserved quan- 
tum numbers, we shall treat all physical transitions related to each separate state n as the 
independent events. 

The exact solutions of a single-particle Dirac equation in external fields can be used for 
calculating the scattering and pair production probabilities in quantum field theory [61 [7J [8]. 
Consider the in- and out-going asymptotic solutions ±i/j(z) for one-particle states given in 
Eq. (jUJ as well as the in- and out-going asymptotic solutions ± 'ip{z) for one-antiparticle states 
given in Eq. (jlBl) . Their asymptotic forms are found after solving the Dirac equation (j5J in the 
potential (jSJ) exactly. As a result, all these in- and out-going one-particle (antiparticle) states 
are related by Eqs. (I50P and (15T1) with the coefficients /(A) and R(X) and their complex conju- 
gated given in Eq. (j49l) . and also in Eqs. (133|) and (JMp . Accordingly, any field operator can be 
expanded either over all possible in-states or, equivalently, over all possible out-states in terms 
of creation and annihilation operators. These are also separated as in- and out-operators in 
such a way that the operator a™/a° ut annihilates the in/out-going particle with the wave func- 
tion + ?p(z) / _?p(z), while the operator (6™) / (6^ ut )^ creates the in/out-going antiparticle with 
the wave function + ip(z)/~ip(z), respectively. The key point [5J [7J [S] is now that the asymp- 
totic relations (150P and (1BT1) between + ip(z) / _ip(z) and + ip(z) (z) determine completely the 
coefficients of the Bogoliubov transformation between the in- and out-going operators 

H^) nut 1 /MUJAt /Wn\t _ 1 nut ^ W 



a n ~ r>(\ \ a n 



(Cf , = ~^TS < Ut - (Cf , (68) 



R(X) n R(X) v n 1 ' v n} R(X) n R(X) 
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whose anticommutation relations preserve the current conservation (1561) for supercritical po- 
tentials. Thus the Klein paradoxical result |.R(A)| 2 > 1 of a single-particle treatment is a 
direct consequence of the fact that fermions obey the Pauli principle and therefore quantized 
by anticommutators [TBI IT7] . 

With these operators, one can construct two complete sets of in- and out-states from the two 
vacua. These are the in- vacuum state |0™) and the out- vacuum state |0° ut ) related for a certain 
n by the unitary ^-matrix operator S n as follows |0Jf) = if> n |0£ ut ). The both are annihilated 
by the operators a™ and a° ut , respectively. This allows one to calculate all necessary matrix 
elements. For example, the absolute probability amplitude for pair creation in a given state is 



w 



pair 



/pout Iflut „out 
\ U n In n 



/(A) 



01 



01 



/(A) 



<0rK n ), (69) 



where (0° ut | 0™) is the probability amplitude for the vacuum to remain a vacuum. The absolute 
probability for production of a single electron-positron pair in a certain state is therefore 



\w 



pair 1 2 



u i (or i on i- 



(70) 



where ti is the relative probability ( 165]) for production of a single electron-positron pair in 
the state i and |(0° ut | 0J n )| 2 is the probability that no pairs are created in this state, or the 
vacuum-to-vacuum probability. This probability can be determined consistently as follows. 

First, we note that the absolute probability ([TO"]) must be equal to the mean number of 
pairs created from the in- vacuum in a given state. Taking the expectation value of the number 
operator, we calculate this explicitly 



fk = ( 



0' 



\R(XW 



0J 



V? (ft*) 1 



ov 



\R(XW 



u 



(71) 



with the result as in Eqs. ( 1661) and ( 167|) . where the transmission coefficient U and the reflection 
coefficient are both referred to a certain state i. From Eqs. ( |70|) and (JTTj) it follows that 



i(orioni 2 = i/r 4 . 



(72) 



Second, we use the probability conservation in a given state together with the Pauli principle 
for fermions to obtain directly, 



i(orioni 2 =r 4 Kori onr 



(73) 



where we apply Eq. ( 159]) for supercritical potentials. This defines the vacuum-to-vacuum prob- 
ability |(0° ut | 0j n )| 2 as in Eq. (1721) . 

Finally, we calculate the absolute probability amplitude for elastic scattering in a given state 
n as follows 



w 



scatt 



out 



„out 



R*(X) 
R*(X) 



K n ) t 

gout 



01 



On 

< Ut « U *) t 
0n) + 



,R*(A) 



w 



pair 



R*(X) 
/(A) 



nout l„out iput I nin\ 
U n \ a n U n \ U n / 



<or |oir> 



(74) 
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where we use Eq. (169]) and the current conservation fl56l) for supercritical potentials. With 
Eq. ( 157]) . the absolute probability for elastic scattering of electron in a certain state is therefore 



where the total reflection of electron by supercritical potentials is due to the Pauli principle 
since the initial state is already occupied. 

We conclude therefore that all Eqs. (172]) . (173]) and ( 175]) define consistently the vacuum-to- 
vacuum probability for a certain z-state as 1/rj, where rj is the reflection coefficient floT]) referred 
to this state. For all these states, the total vacuum-to-vacuum probability is 

|(0 out |0 in )|2 = "Q (ri) -i = exp f-^mr.j , (76) 

where the products and the sums are taken over all states i = (p, a) with the quantum numbers 
p = {e, p_i_} and also with the two spin projections a = 1, 2. Since |(0 out |0 in )| 2 < 1, the vacuum 
is unstable producing particle-antiparticle pairs under the influence of supercritical potentials. 

The total probability for vacuum decay via pair creation reads 



W = 1 - exp I - V" In n J ~ ln r i = 2 ln Ti ' ( 77 ) 

\ i / i p 

where the factor 2 is due to summation over two spin projections a = 1,2. The remaining sum 
over p = {e, p_i_} can be rewritten as the integral in a box-like volume V±T with V± — J d 2 x± = 
j dx\dx2 being the area transverse to the ^-direction, and T the total time. This yields 

where r = r(e, p±) with p± = \p±\ being the rotationally invariant function as it follows from 
Eq. ( IBT]) . This invariance reduces the integral J d 2 p± to n J dp\. For supercritical potentials 
with v > m, integrating over e and p\ has to be done over the Klein region 



p ( ^ = e + v > ^p\ + m 2 , p^ = e — v < —^p 2 ± +m 2 . (79) 

The total probability is now completely defined by Eqs. ( 177]) - (179l . The vacuum decay 
rate (the total probability per unit cross- sectional area and unit time) is represented due to 
Nikishov [7J as 



w± = fo , 2h3 / dp 2 ± / de\nr(e,p 2 ± ), (80) 

where the integration region in the (p^_, e)-plane is shown in Fig. 1. The reflection coefficient 
r { £ iP'±) i s defined by the symmetric under the interchanging /j, •<->■ v expression ( |57l) . and 
therefore remains invariant upon replacing which interchanges merely fi and v given 
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v — m 




—v + m 

Figure 1: In the (jp\, e)-plane, the integration covers the positive region restricted by two 
intersecting parabolas e = v — \/p\ + rn 2 and e = —v + \Jp\ + m 2 with horizontal axes of 
symmetry above and below the p^-axis for v > m. 



by Eqs. ffl2|) and (TT9T) . This allows us to obtain from flHUl) yet another representation for the 
vacuum decay rate as follows. We interchange the order of integration 



1 



w± 



(2vr) 2 ft 3 



Me+v) 2 -rn 2 

de l dp\ 

—v+m Jo 



v—m 

de 

o </o 



(e— v) 2 — m 2 



dpi 



In r(e,p 2 ± ) 



(81) 



and replace e — > — £ in the first term of Eq. (1811) . Due to the symmetry r(— = r(e,pj_), 
the two terms in Eq. f|8T|) are then combined into a simple integral representation 



1 



2n 2 h 3 



de 



(e—v)—m 



dp\ In r(e,p 2 ± ) 



(82) 



where the integraion region in the (e,p5_)-plane is shown in Fig. 2. In the same way, we express 
also the total mean number of created electron-positron pairs 



n 



v ■ 



(83) 



where the average number of electron-positron pairs created in a certain z-state is defined 
by Eq. (166|) independently on the spin state a = 1,2. Like the reflection probability, it is the 
rotationally invariant function n = n(e, p±) with the same symmetry n(— £,p^_) = n(e,p 2 ± ) 
under the interchanging \i -<->• v in Eq. ( |66l) . The corresponding integral representation is 




Figure 2: In the (e, p^_)-plane, the integration covers the region under the left branch of the 
parabola p\ 



le-v) 



m in the first quadrant for v > m. 
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therefore obtained as before via replacing the sum over p = {e, p_i_} in Eqs. f l83|) by the integral 
over the Klein region. This yields the mean number of electron-positron pairs produced by 
the supercritical Sauter potential (J2]) from the vacuum in unit area per unit time, i.e., the pair 
production rate, 

n± = J de J dp\ n(e, p\) , (84) 

with the integration region in the (e,Pj_)-plane shown in Fig. 2. 

Eqs. (I82p and ( 154"]) allow us to calculate the pair production rates for supercritical poten- 
tials with v > m. In the limit k — > and v — > oo (the linear potential due to constant 
electric field), the obtained results can be compared with the Schwinger formula (JTJ, and in 
the regime close to this limit, with semiclassical expressions of Refs. [15] H6j S3 SHI SHI 150] . 
The opposite limit k — > oo (the step potential due to delta-shaped electric field) is however 
beyond the semiclassical approximation. It leads to infinite results caused by the extremely 
strong electric field at the origin [17] . In the regime close to this limit, no analytic results for 
production rate were obtained so far to the best of our knowledge. 

In equation fl82|) . the reflection coefficient is given by Eq. (1ST]) , where /i = p,(e,p±) and 
v = u(e,p 2 ± ) are the functions of e and p\ defined by Eqs. ffl9l) . ffT2|) and pip with the 
constraint ()2Up . while A = v / kh > is a constant. Taking logarithms of Eq. (|57p leads to the 
expansion 

lnr = 4^^ — cosh(27mA) sinh(27rn/i) sinh(27rnz/) , (85) 
n=i n 

where the right hand side is found by replacing each logarithm containing the ratio of hyperbolic 
functions as ln(sinhx/ sinhy) = — J2^=i cosh(2nx)/n + Y^n=i cosh(2ny)/n, and combining all 
sums. 

With Eq. ( 185|) . the vacuum decay rate (1521 takes the form 

2 



w± 



.,.„ Y - cosh(27rnA) , (86) 

n=l 

where are the integrals 



pv—m rv—m f(e—v) 2 —m 2 

J (n) = / del {n \e)= del dp\ sinh2vrn/i(e,^) sinh27rnz/(e,^) , (87) 

Jo Jo Jo 

with the integration region shown in Fig. 2. 

The ^-integration in Eq. ( 157)) can be done exactly [53]. The result is expressed in terms of 
the total energy e with the help of the functions 

0±(e) = /i(e,O)±i/(e,O), 0+(e) 0_(e) = (ve)/(hk) 2 , (88) 

where /i(e, 0) and u(e, 0) are two positive square roots 

H{e, 0) = y/(e + v) 2 -m 2 /2hk , i/(e, 0) = ^{e - v) 2 - m 2 /2hk , (89) 
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as in Eqs. f lT9|) . ffl~2]) and ffTTj) but now with p 2 = 0. Using these functions, we obtain 

= F^(6+(e)) - F {n \64e)) , (90) 



where 



Fto { 0) = ( "Vf(»V)-M!^Vf<") (9 ), (91) 



27m/ 1 w 2 V 



with 



F[ n) (9) = 2irn 9 sinh (27m 9) - cosh (27m 9) , (92) 
= Chi (27m g) - 2nn6 Sinh (W) + 2 COSh (27r " g) , (93) 

where Chi (2irn9) is the hyperbolic cosine integral and the next two terms represent the leading 
terms of its asymptotic expansion for large arguments 27m 9, n = 1, 2, . . . . 



Now the remaining e-integral in Eq. ( la?) reads 

pv—m 

J^= / de&\e) = jj n) - J W , (94) 
Jo 

where 

rv—m 

4 n) = / deF^ n \9 ± {e)). (95) 



The two terms in Eq. ( 19~4"|) can be combined into a single integral by subjecting J±^ in Eq. (195]) 
to the following change of variable 

f m 2 1 1/2 

£ W= M T" [^-(^p] ) • (96) 



with < 9 < \Jv 2 — m 2 /hk, where the endpoints are the zeros of the function e(9), whereas its 
maximum (v — m) lies at y/v(v — m)/hk. Note that m 2 < [v 2 — (hk9) 2 } within these limits. 



For < 9 < yv{v — m)/hk, where the function ( |96|) increases from to (v — m), the positive 
square roots ( 189]) are expressed in terms of 9 as follows 

V(£±vr -- m ^ ± M0 + J {1 - [v2 ^l ker] ) 1,2 >0 (97) 

and therefore, 

9-{e) = 9, 9 + {e) =ve{9)/{hk) 2 9. (98) 

Within this range, we substitute 9 instead of e into the second integral in Eq. ( 1951) . This 
yields 



J W =/ deF^\9^e))= d 9^-F^(9). (99) 
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In turn, for \Jv(v — mj/hk < 6 < y/v 2 — m 2 /hk, where the function (1961) decreases from 
(v — m) to 0, the two manifestly positive square roots (1891) are 



f m 2 ^ 1/2 



(100) 



Combining these yields 

e + (e) = 6, 6 -(e) =ve{8)/{hk) 2 6 . (101) 
Within this range, 6 replaces e in the first integral in Eq. fl95|) . This gives 

jj n) = / deF^\6+(e)) = - / dO^-F<- n \0) . (102) 

With Eqs. (11021) and ( 1991) we obtain the integral ( 194"1) in terms of dimensionless variable 9 as 
follows 



Jo dd J 



27rn/ d# 6 \ hk J dd 



, (103) 



where 



= Vv 2 -m 2 /hk 1 (104) 



and the last term is found via partial integration thanks to the vanishing of the function fl96l) 
at the endpoints. Substituting (152]) and (l93l) into Eq. (I103f) yields 

r S 

jW = (ft*;) 3 / d9f(6) cosh(2nn9) . (105) 

Here the function f(9) has the form 

f(6) = (Hk)- 2 [(hk6) s(6) - (v 2 /3) e\6)/(hkef] , (106) 
and by Eq. (|94"|) reads explicitly, 

(/)2 /J2 \ 1 / 2 \2 / n2 o2\ 3 / 2 
-H^) <io7) 

Note that the integral over this function vanishes j^dO f(6) = 0, so that J(°) = 0. The integrals 
are all functions of v,m, and k. 

With integrals f ll05[) the vacuum decay rate (1861) becomes 

2k 3 f s ^ 1 

w ± = — d6f(6) > - cosh(2vrnA) cosh(2vrn0) , (108) 

where we interchange the order of summation and integration. For < 9 < 6 < A, the sum 

CO . ^ CO ^ 

^ - cosh(27mA) cosh(27rn0) = - - (e 2nnX + e ~ 2 ™ x ) cosh(27rn0) . (109) 



n 2 ^— ' n 

n=l n=l 
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can be found as follows. The second sum in Eq. (I109j) has the form 

™p n l 

y — cosh(nx) = — - ln(l — 2pcoshx + p 2 ) , p 2 < 1 , (HO) 

n=l 

with x = 2tt9 and p = exp(— 2n\). Explicitly, it reads 

e" 2 ™ A cosh(27m0) = -i In [l - e - 2 -( A + e )] - i l n [l - e -MA-*)] . 

n=l 

The first sum in Eq. (11091) can formally be written as 

- J] - e 27mA cosh(27m0) = ± Li^e 27 ^) + ± Li^e 27 ^) , (112) 

n=l 

where Lii(z) is the polylogarithm function 

00 n 

L^)=Ei> 1^1 <!> ( 113 ) 



n=l 



with z/ = 1 for which we use the analytic continuation into the region \z\ > 1, 

U 1 (z) = U 1 (-) -ln(-z) . (114) 



Applying this to Eq. ( 11121) . we obtain 

^ E cosh(27m£) = -ttA - ± ln [l - e ~ 2 < x ^] - - ln [l - e" 2 ^-*)] . (115) 

n=l 

By Eqs. (II 111) and (11151) . the infinite sum (1 109f) becomes 

00 

J3 - cosh(2vmA) cosh(27r^) = -ttA - - ln [l - e - 2 ^ x+9) ] - - ln [l - e - 2 -( A - e )] . (116) 

n=l 71 

Finally, inserting Eq. (11161) into Eq. (I108D yields the vacuum decay rate 

~ Jj0f(9)\n[l-e- 2 ^}, (117) 



w± 



where we use the symmetry of the function f{9) = /(—#). 

The pair production rate (1541 can be obtained in the same way for large A, i.e., close to the 
constant-field limit. In this case, the average number of produced pairs (1661) becomes 

2 sinh 2ttu sinh 2ttv , n . , , 

n = -. — —- -, — ~ 4 e~ 2nX sinh 2tt/x sinh 2txv . (118) 

[cosh 2-7rA — cosh 27r(/x — v)\ 

Substituting Eq. (jTTgjl into Eq. (JHJD yields 

_ 2 _o„\ 2k^ 



n± 



n 2 K 3, 



e ~2nx j(i) = e - 2 .A / d 9 f{9) cosh{2ir9) = ^ / d9 f [B)e 2 < e - X \ (119) 
71 Jo 71 J-e 
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where the integral J^ 1 ' is determined by Eqs. ( 18 7p and f 1 1 U 5 j) for n = 1. The obtained equa- 
tions (I117p and (11191) provide us with simple integral representations for probabilities ([77]) 
and ( 183|) . respectively, where the sum over all quantum numbers a, po and p_i_ is replaced by 
the ^-integral with the function f(6). In particular, the pair production rate fl 11 9 [) is given by 
the first term in the expansion of logarithm in Eq. (JTTTJ) for the vacuum decay rate, i.e., in the 
same way as Eq. ( 153"]) is related to Eq. ( 177|) . 



4 Pair production rates and field inhomogeneity 

The obtained rates ( JTT7]) and (1119j) include the effect of spatial variations of the electric field (J2J) 
on pair production. To make this more explicit, we represent the function ( I107p in the form 
f(0) = (-l/3)dg(0)/d0 with 

9(0) = o ( \ 2 _ e J )l/2 = -g(-e) , (120) 

and perform the partial integration in Eqs. ( 11171) and (11191) by taking into account that the 
function (I120p vanishes on both ends. Then the vacuum decay rate (11171) becomes 

2k 3 r s 1 

»-¥/*' |, '^ry (121) 

In what follows, we shall use the natural units in which the height v and the width 1 jk of the 
potential barrier ([3]) are measured in units of the rest energy mc 2 and the Compton wavelength 
X e of electron, respectively 

v = amc 2 , l/k = (3X e , (122) 

where a and (3 are dimensionless parameters, whose ratio is 

e = a/ (3 = E /E c . (123) 

Note that a > 1 and therefore /3 > 1/e for supercritical potentials. By these restrictions the 
electric field (T5]), whose peak E is of the order of presently available field strengths e <C 1 
can only create the pairs near the constant limit /} > 1, whereas for pair production near the 
sharp limit /3 <C 1 much stronger field strengths e ^> 1 are required. In the units ( 11221) . the 
parameters (121!) and (I104p entered in Eq. (I12ip read 

\ = a/3, 6 = f3Va 2 - 1 . (124) 

Setting these up in Eq. ( 11211) and expressing the integral in terms of the new dimensionless 
variable $ = a — 9 / (3 yields the rate per area 

We define now the vacuum decay rate per volume as w ee kw± = w±/ (3X e . It can be 
obtained by dividing both sides of Eq. (I125p by the width 1/k of the electric field §2§. This 
yields explicitly, 
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Figure 3: The dimensionless rate (11261) is plotted as a function of two independent variables in 
natural units: the width (3 and the height a. 



The expression fl 1 2 6 [) is a function of two independent variables. Similarly to the dimensionless 
field (I123p . it increases with growing a and decreases with growing (3 as shown in Fig. 3. For a 
fixed, the rate H126[) vanishes in the limit (3 — > oo and diverges in the limit /3 — > 0. In the later 
limit, the electric field takes a singular 5-shaped form at the origin. The divergence of the 
vacuum decay rate (11261) in this limit is the property particular to the sharp step [17] . whose 
discontinuity leads to enormous electric field. 

With increasing both arguments to infinity in such a way that their ratio e remains a constant, 
the expression (11261) becomes a homogeneous function of this ratio. It is therefore also a constant 
in this limit which we denote by w M (e). The rate u> lcf (e) corresponds to a spatially uniform 
electric field eE c obtained from the Sauter field ([2]) in the same limit. For its explicit calculation, 
we rewrite Eq. (11261) in the form 

w{a > « = £s <*C^»( i - !) (*-£)"*(*- 1 - 9 V^i 

where w = 2Ti(3 r d. Substituting then a = e(3 into Eq. (I127P and taking the limit ft — > oo with e 
fixed, we obtain 

w M (e) = w(a = eP,P)\^ 00 = 6 2 / d,," 1 ' 2 ( W - vr/ e ) 3 / 2 — — . (128) 

The integral (11281) can be evaluated explicitly by substituting into Eq. (11 28[) the expansion 




(129) 
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(a) (b) (c) 

Figure 4: The dimensionless rate (I126P is plotted as a function of the field width in Compton 
units for fixed values of the dimensionless electric field e = 0.1 (a), 1.0 (b) and 10.0 (c) (solid 
blue lines). Dashed red lines refer to the constant-field rate (11301) for these values of the field. 
The variable /3 starts running from 1/e. 

and interchanging the order of summation and integration. In this way, we obtain the constant- 
field rate per volume as 



w {l > = wxT 

c n=l 



CO f*C0 

1 J ir/e 



7r/e 



c ^ e-"^ fl tvk_\ 
47r 5/2 A 4 e 2^ n 2 u I 2 ' ' e J 



where [7(1/2, — 1, z) is the Tricomi confluent hypergeometric function [58] with the argument 
z = rwr/e and the first term represents the constant-field limit of the total mean number 
of created pairs (I119p . Thus, the constant-field limit of the vacuum decay rate (I126p is the 
expression (I130p rather than the Schwinger formula ([1]). It agrees, however, with formula (JT]), 
if the latter is extended to the space-dependent electric field (J5J) and is averaged over the 
infinite width of a spatial variation appropriated for a constant field. The expression (11301) 
represents therefore the locally constant-field rate. We explain this fact and present the detailed 
comparison of the two constant-field rates in Appendix A. In formula (I130p . the ratio e takes 
on arbitrary values despite of the conditions /J > 1 and a ^> 1. The vacuum decay rate 



approaches the constant-field limit (I130p with increasing j3 and fixed e as shown in Fig. 4. 

Below the constant-field limit (I130p . the expression (I126p for the vacuum decay rate per 
volume interpolates analytically between the regime of sharp field /3 < 1 with e 3> 1 and the 
regime of constant field /3 > 1 with an arbitrary e > 1/(3. In order to describe such a behavior, 
we put the integral (I126P in a more symmetric form by noting that the upper and lower limits 
of integration are the two zeros of the quadratic polynomial in the numerator. By denoting 
these limits as r d± = a ± y/a 2 — 1 with $ + -d- = 1 and + = 2a, we represent the rate (11261) 
as follows 

w(a, (3) = f + dd [(* + +4_)/2-0\ [(^+-^-^-)] 3/2 (131) 

3vrAfV J p + -tf)(tf-tf_)+TML] 1/2 e^-l' 

where d + > 1 and < 1. The integral (I13ip is dominated by the region $ < 1/27T/3 in 
which the function 1/ [e 217 ^ — l) differs appreciably from zero. With respect to this region, 
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Figure 5: The interval of integration with endpoints $± = a ± \/ a 2 — 1 is shown with respect 
to dominant region of the width 1/(3 in Eq. f 1 1 3 lj) : completely outside (a), partly outside and 
partly inside (b), completely inside (c). The corresponding analytic expressions are given by 
Eqs. ffT35l) . f TltoD and ffT52l) . respectively 



the positions of endpoints $± are completely arbitrary. Thus locating the interval of integration 
as shown in Fig. 5 we can estimate the integral (I13ip analytically as follows. 

We assume first that the position of the upper limit is far to the right of the dominant 
region i? + 3> 1/(3 as shown in Figs. 5a and 5b. Then the region of integration is cut off by the 
factor 1/ (e 2 ""^ — l) and is of the order 1/(3. In this region, i? <C $ + and the upper limit can 
be extended to the positive infinity. Thus the rate (11311) becomes approximately, 



00+ >1 d7rAe ^- e p -1 

POD 1 

A- 2 / duo u- 1/2 (w - A) 3/2 = iu(A) , (132) 



3ttA 4 



where we replace the variable $ by w = 2n(3d and introduce the scaling parameter A = 27r 
The integral is now performed by substituting dl29[) into HI 32 j) and interchanging the order of 
summation and integration. Then we obtain 



n=l 

c ,_.^e~ nA 



Air 1 / 2 ** 



^E^T^^A-UA) , (133) 



71=1 



where all steps of calculation are similar to those in Eqs. (I128l) - (ll30p . This leads to the same 
hypergeometric function [7(1/2, — l,z) but now with the argument z = nA. As a such, the 
parameter A depends on a and (3. It can also be expressed in terms of a and e with the help 
of Eq. ffT23]) . This yields A = (vr/e) (2a/tf+) with A > (vr/e). Thus, the rate ffT33D is always 
smaller than its constant-filed limit (I13(jp . The latter is attained by Eq. f 1 1 3 3 j) in the limit 
(3 — >■ oo with fixed e via substituting A — > ir/e. 

In going from Eq. (I13ip to Eq. (11331) the position of the lower limit of integration is still 
arbitrary. For its location there are two possibilities: far to the right 3> 1/(3 (Fig. 5a) and 



23 



far to the left <C 1/(3 (Fig. 5b) of the boundary 1//3. For the first, the interval of integration 
is located completely outside the dominant region, whereas for the second, it is partly outside 
and mostly inside. The both are the limiting cases of Eq. f 11 3 3 j) and can be estimated as follows. 

If the interval of integration is located completely outside the dominant region (Fig. 5a), 
the rate (I133p is obviously very small. The corresponding restriction 3> or (3d- 3> 1 
leads to (3 ^> d + > a > 1 > This implies, in turn, that the initially imposed restriction 
$ + ^> 1/(3 is already preserved, and also that the electric field ([2]) is near its constant limit 
1 with |3 > a or e < 1. As long as A ^> 1, we estimate the rate fl 1 3 3 [) by making use of 
the asymptotic expansion of the function [7(1/2, —1, z) for large arguments z = nA [58] : 

[7(1/2, -1, z) ~ z~ 1/2 [l - (5/A)z- 1 + O (z~ 2 )] , z ->• oo . (134) 

Inserting this into Eq. HI 33j) and subjecting further to A ^> 1 yields 



w(A) 



A > 1 



47TV2A4 



A -5/2 



Lis (e 
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- A ^-^A- 1 Li 7 (e' A ) + 
4 2 v 7 



, 4vr 1 / 2 A 4 
A > 1 e 



A -5/2 



cxp 



-A( 1 + - A~ 2 + (A~ 3 ) 



(135) 



The rate f 1 1 3 5 f) describes the pair production in the weak-field regime 6 C 1 near the constant 
limit ^> 1 with 1 < a ^ 0. In terms of a and e, it can be represented via the substitution 

A = (vr/e) (2c0?+) as 



muck, e 



c 5/2 



^ / Q 47T 3 Af 

e < «/»+ e 



\ 5/2 I 7r2a 

2^ J exp ri^ 
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7r 2a 



(136) 



where e <C «/■$+ < 1. The leading term of the expansion (I136P with slightly different pre- 
exponential factor including \J a 2 — 1/a instead of $ + /2a was found in Ref. [H] by integrating 
the Nikishov result and also in Refs. [15J H7J HH [50] by various semiclassical approximations. □ 
The corrections in powers of small e were obtained in Ref. |49j . 

In the limit — > 00 with fixed eClwe obtain from (I135P the small-e asymptotic (I17ip of the 
constant-field rate f )130p by replacing A — > 7r/e 3> 1- With increasing a from moderate a > 1 
to large a>l values until a <C (3, the rate fll35p increases approaching the asymptotic f )17ip 
from below. For a 3 3> 0, the corrections can be obtained by expanding the scaling parameter 
A ~ (tt/c) (1 + l/4a 2 + • ■ ■) and then the expression ( I135P as follows 



47T 3 )l 4 

e < 1 47F Ae 



e 5/2 exp 



7T 



5e 2 1 ^ 
4 7r 2 4a 2 



e 2 a 4 



a < /3 < a 3 



(137) 



If the interval of integration is located only partly outside the dominant region while oc- 
cupying the most of its part inside as shown in Fig. 5b, the rate (I133P becomes very large. 



1 Note that our parameter a = v/mc 2 > 1 coincides with I/7 of [45j [47], with 1/e of [48l [49] and with a 
of [SO]. 
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The corresponding restriction <C 1//3 together with the initial restriction $ + ^> leads 
to the condition a 3> 1 for which the two restrictions are reduced to a/3 3> 1 and (3 /a <C 1, 
or e > 1. As long as A C 1, estimating the rate (I133p becomes more involved than in the 
previous case. Indeed, the use of the asymptotic expansion of the function [7(1/2, —1, z) for 
small arguments [58] : 

[7(1/2, -1, z) ~ ^= (l - | + O (z 2 )) , (138) 

with z = nA leads to a slow convergent series for large n in Eq. (11331) . Therefore we determine 
the small-A behavior of the rate (11331) from the integral representation (I132p as follows. First 
of all, we rewrite this in the form 

c ,_ a r , a aV /2 co 



w(A) = — — A" / dull--] . (139) 

1 ' SvrAf A V u) e»-l y } 



In order to perform the integral (11391) for A <C 1, we expand 

u) 



A = f; a n A n U- n = 1 - ^uT 1 + ^V 2 + ■ ■ ■ , (140) 



n=0 

with the coefficients 

«n = — ; — w — r-r • (141) 

2 n (2n- l)(2n-3)n! 

Substituting (I140p into (11391) and separating the integral into two parts yields 

„(A) = S- A- V a „A» (r^l- / A . (142) 

Here the first integral is the product r(2 — n)C(2 — n), where r(2 — n) and C(2— n) are the gamma 
and zeta functions, respectively [SB]- In the presence of T(2 — n), each term of the infinite sum 
with n > 2 over this product consists of the singular part (— l) n ((2 — n) /(n—2)\ = B n ^\j (n— 1)! 
plus the regular part (— l) n tp(n — 1)C(2 — n)/(n — 2)! = t/>(ti — l)S n _i/ (n — 1)!, where -B„_i are 
the Bernoulli numbers and ip(n — 1) is the polygamma function [58]. In the second integral of 
Eq. (11421) with the small upper limit A < 1, we substitute the expansion 

— = Eiif'-' = 1 -2 + i2--' (143) 

k=0 

and extract from the double sum the logarithmically divergent term with n = k + 1 involving 
the coefficients B n _i/(n — 1)! after the integration. This cancels the singular part coming from 
the gamma function T(2 — n) to each power of small A. Then, the rate (11421) takes the form 



w(A) = ^rrA- 2 
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(144) 
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For A 1, the expansion (I144p is well approximated by the first few terms. In doing so, we 
obtain finally, 



iu(A) 



3vrAf 



■A" 



A < 1 
where 7 ~ 0.5772157 

Si = 

So = 
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-lnA + ^x 



+ A 2 



S 9 



+ O (A 3 ) } , (145) 



is the Euler constant [58] and the sums Si and S2 read explicitly, 



1 OO 



n=0 
n^2 
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1 OO 



n 



2) 



0.5411802 



(146) 



The rate (j!47p describes the pair production by extremely strong electric fields e ^> 1 within 
the following range of the parameters a ^> 1 and 1/a <C /J < a. Under these conditions, the 
electric field (j2J) can either be near the sharp limit 1/a < /3 < 1, or near the constant limit 
1 < ^ < a. Thus, the two regimes of pair production ft >C 1 and /5 > 1 bear a neat similarity 
for very large values a ^> 1. 

In the limit ft 00 with fixed e > 1, we obtain from (I145P the large-e asymptotic (11771) of 
the constant-field rate (I130p via substituting A — > 7r/e <C 1. As long as a > 1, the parameter 
A in (I145P can be expanded as A ~ (tt/c) (1 + l/4a 2 + • • •) and the expression (I145j) can be 
then re-expanded in powers of small 1/a as follows 



w(a, ft) 
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l/a</3<a. (147) 



In the regime 1/a < ^ < 1, the rate (I147p is very large but is overestimated by the large-e 
asymptotic fl 177ft . With growing ft the expansion (I147p decreases approaching the asymp- 
totic (I177P in the regime 1 <C /5 < a. 

Let us assume now that the position of the upper limit lies far to the left $ + <C 1/ft 
of the boundary 1/ft as shown in Fig. 5c. The restriction i? + <C 1/ft or ftd + <C 1 leads to 
ft «C 1?- < 1 < a < $ + . This yields, in turn, the restriction <C 1 for the lower limit 
and implies that the electric field fl2]) is near the sharp limit ft <C 1 with ^ < a or e > 1. 
The parameter a is restricted by the condition « < < 1//?. The integral (I13ip . where the 
interval of integration is located completely inside the dominant region can be estimated as 
follows. Within this interval, we use the representation 

2-nftd 

k\ 



1 



OO R 

(2nft$) k - 



k=0 



2nft$ 
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where are the Bernoulli numbers 
the change of integration variable d = 



Substituting Eq. (I148P into Eq. (I13ip and making 
- v + (t? + + #_)] /2, we obtain 



w(a, /3) 



3vr 2 Af /3 
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(149) 
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where a = \/a 2 — 1/a < 1 and the integrals 



( (o) = a'Jjv (1 V _ a2 J )1/2 (1 - H fc 



(150) 



are certain combinations of the complete elliptic integrals K(o) and E(a) of the first and the 
second kind, respectively. Explicitly, these read 

J (a) = - [(8 - 7a 2 ) E(a) - (8 - 11a 2 + 3a 4 ) K(a)] , h(a) = , 
3 

J 2 (a) = — [(8 + 5a 2 + 3a 4 ) E(a) - (8 + 9a 2 + 12a 4 ) K(a)] , . . . . (151) 

For a/3 -C 1, retaining only the first few terms in Eq. (1149 j) provides us with a good approximate 
rate. In fact, the expression (11491) represents a quickly convergent series 



w(a, (3) 



a 3 



(3d + < 1 37r ' A ' 13 
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I (a) + -/ 2 (a) (a/3) 2 + ((a/3) 3 ) 



a/3 < 1 . (152) 



The rate (I152p describes the pair production in the strong-field regime e ^> 1 near the sharp- 
field limit /3 <C 1 for all values of the parameter a satisfying 1 < a <C 1//3. With increasing a 
from moderate a > 1 to large a ^> 1 values until a 1//3, it can be expanded further as 



2c a 3 

/3 < 1 J7r A e P 



l<a<l//3, (153) 



where a 3 //3 = ea 2 <C e 2 as long as a/3 <C 1. The rate (I153p is therefore smaller than the 
rate (I147p in the same strong-field regime of pair production. 

The pair production near the constant-field limit (3 3> 1 is now described as follows. For 
(3 ^> 1 > the restriction (3d + 3> 1 is always satisfied. For moderate values of the parameter 
a > 1, we obtain in addition the restriction /3i9_ ^> 1 due to (3 ^> 1. This yields (3 ^> 
?9 + > a > 1 implying that the regime of pair production is just the weak-field regime e < 1 
described by the rate (I135p . With growing a from moderate a > 1 to large a ^> 1 values 
until a <C (3, the first restriction /3$ + ~ a/3 ^> 1 is satisfied automatically, whereas the second 
/3$_ ~ (3/ a ^> 1 is still preserved by the condition 1 < a < /3. With these conditions, the 
weak-field regime is still maintained but is described now by the large-a asymptotic (11371) of 
the expansion (11351) approaching from below the small-e asymptotic (I17ip of the constant-field 
rate (11301) with growing a. Increasing a further as a ^> (3 ^> 1 leads to /3$_ ~ /3/a 1 instead 
of (3i9_ ~ (3 /a ^> 1. Then the pair production undergoes the transition from weak-field e <C 1 
to strong-field e > 1 regime, where it described by the large-a expansion (11471) approaching 
the large-e asymptotic (I177P of the constant-field rate (11301) . The increase of pair production 
rate due to transition is of the order of two magnitudes. Near e = 1, the rate of pair production 
is interpolated analytically by the expression (I13ip . 

Consider now the pair production near the sharp- field limit /3 < 1. For /3 <C 1 < the 
restriction [3d- <C 1 is always satisfied. For moderate values of the parameter a > 1, there 
is in addition the restriction /3i9 + <C 1 due to (3 <C 1. This yields (3 <C < 1 < a. The 
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regime of pair production is therefore the strong-field regime e ^> 1 described by the rate (I152p . 
With increasing a from moderate a > 1 to large a ^> 1 values until a <C 1//3, the first 
restriction ~ (3 /a <C 1 is satisfied automatically, whereas the second (3$ + ~ a (3 <C 1 is 
still preserved by the condition 1 C a C The strong-field regime e ^> 1 of pair production 
in now described by the large-a expansion (I153p . Increasing a further as a > > 1 yields 
/3$+ ~ a/) > 1 instead of (3{} + ~ aj5 <C 1. Then the strong-field regime of pair production 
undergoes the smooth transition from e 1//3 2 to e ^> 1//3 2 , where it described by the large-a 
expansion (11471) . 

For very large values «>1, the pair production is thus described by the expression (11471) 
whatever ^> 1 or /3 <C 1. In the strong-field regime e> 1, this expression becomes very large 
but does not exceed the large-e asymptotic (11771) of the locally constant-field rate (11301) . The 
latter is therefore the upper limit for the vacuum decay rate (11311) . Thus, in agreement with 
Refs. [331 113 S3 EH El], we come to the conclusion that the spatial variations of the Sauter 
field (T5]) cannot increase the pair production over the locally constant-field rate (11301) even 
for extremely strong fields with e ^> 1. The later rate overestimates however the Schwinger 
formula f !159j) as long as e 3> 1. 

5 Conclusion 

We have calculated analytic expressions for the production rate of electron-positron pairs from 
the vacuum by the Sauter potential. For an arbitrary potential barrier, the rate was related 
to the scattering amplitude on the barrier, and expressed as an energy-momentum integral 
over the logarithm of the reflection coefficient. For the Sauter potential, we have evaluated 
the rotationally invariant integral over transverse momenta in three dimensions exactly. The 
remaining integral over the energy gave us the simple spectral formulas for the vacuum decay 
and pair production rate. The analytic expressions for these rates were derived for the entire 
range of the parameters v and k of the supercritical potential ()3]). This allowed us to access 
different physical regimes from weak to strong fields. 
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In previous studies of pair creation by inhomogeneous fields, the transition to the constant 
field has often been made for the local probabilities before the energy-momentum integration. 
In order to outline this procedure, we consider the densities for reflection and transmission 
probabilities, and also for the average number of pairs created in each mode given by Eqs. f loTj) . 
( 165]) and ( )66|) . respectively. 

The reflection and transmission coefficients in Eqs. ( 15"T|) and ( 16"5"j) can be written as 



A Appendix: Constant-field limit 
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with r — t = 1 for supercritical potentials, where n is the density number of produced pairs ( 166|) . 
In order to find this quantity in the constant-field limit k — >■ with fixed Eq, we first substitute 
v = (\e\E )/k. Then evaluating the limit yields 



n 



exp [-7T (pi + m 2 ) /(ft|e|Eb)] . (155) 



By Eq. (|155p . the reflection coefficient in Eq. (I154p becomes independent of the energy e. Its 
logarithm is 

X 

lnr(p^) = — In (1 — n) = — exp [— rar (p 2 ± + m 2 ) /(h\e\E )~\ . (156) 

n=i n 

With Eq. ( I156p the vacuum decay rate per area (1821) in the constant-field limit becomes 

- ^Es/r <k r* 4e ^ i4,, ' )/ww - <i57) 

Here the p^-integral results in a factor of | e | / ri7r an d the ^-integral, via the substitution 
de = \e\E Q dz, in a factor of \e\E L, where L is the infinite length. This yields the constant-field 
rate per area 

wf(E ) = Lw cl (E Q ), (158) 

where w cf (E ) is the vacuum decay rate per volume given by the Schwinger formula ([!]). For 
the constant electric field Eq, it reads 



n=l e 71=1 



with e = Eq/E c . Thus, the finite rate per volume (11591) is obtained from the infinite rate per 
area (11581) via dividing by the infinite factor of a length. 

In order to compute the limit k — > after the energy-momentum integration, we modify 
the above derivation by noting that E(z) — > Eq for k — > 0, where E(z) is the electric field ([2]). 
Then Eq. (I156p can be rewritten as 

V 

In r (p 2 ± ) = lim'S^ — exp [—nn (p 2 ± + m 2 ) / (h\e\E(z))~\ . (160) 

71=1 

The constant-field rate per area (I157P becomes 
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where we interchange the limit and the energy-momentum integral and perform the integration 
over p\. In the remaining integral, the energy e is related to the position z of the electric field 
E(z) as de = \e\E(z) dz. For the Sauter electric field (j2J), this substitution yields 



^2 e -nn(E c /E(z)) 
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By changing the variable y = kz the constant-field rate per area A162H takes the form 

wf(E )=(hm^jw M (E ), (163) 
where lim(l/fc) plays now the role of the infinite length. The constant-field rate per volume 

k— >0 

reads 
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where x = sinh y. The last integral in Eq. ( 11641) is a certain combination of the modified Bessel 
functions K v {z) with v — 0, 1 and z = n7r/e which can be expressed in terms of the function 
£7(1/2, —1, z) [5H]. This yields, finally, the constant-field rate per volume as 
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This result was already found in Eq. f II 3 H . The rate ( 11650 as well as the rate ( 11591) corresponds 
to the constant electric field Eq but does not coincide with the latter. Moreover, the above 
derivation shows that the Schwinger expression fl 1 5 9 [) cannot be recovered as a constant-field 
limit of the non-constant field rate in which the energy-momentum integral was already per- 
formed. On the other hand, both rates per volume are found from the same (infinite) rate per 
area by the only different factorization of an infinite length and can therefore be related to each 
other. Comparing fl 158[) with ( 11630 . we obtain 

wf(E ) = ( lim \ j w M (E ) = Lw ci (E ) = L lim w ci (E(z)) 

/+oo / i \ r+oo 

dzw ci (E(z)) = (lim- J / dyw cf (E(y/k)). (166) 
\k-+okJ 

This expresses the rate ( 11650 in terms of the rate (11590 extended to the space-dependent electric 
field ()2]) and averaged over the infinite width of a spatial variation 

/+oo r+oo 
dyw ci (E(y/k)) = \imk dzw cf (E(z)). (167) 
J-oc 
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The rate ( 1165)) is therefore the locally constant-field rate. The locally constant field approxi- 
mation was introduced in Ref. [17]. In section 4, we have shown that the expression H165|) is a 
constant-field limit of the vacuum decay rate (11261) . 

The rate (11651) possesses the integral representation (11281) . Similar integral representation 
can also be found for the rate (I159p provided that we express the sum in Eq. ( 1 1 5 9 f) as 
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Inserting this into Eq. (11591) yields 



For better comparison of the two constant-field rates (I159P and (I165p . let us estimate their 
asymptotic behavior for small and large e, respectively. For weak electric fields e<l, a good 
approximation of the rate (11591) is the Heisenberg-Euler expression given by the first term of 
the Schwinger formula 

w«(e) = ^/^(e~^) - ^e'expC-Tr/e) , e« 1. (170) 

The corresponding behavior of the rate (11651) can be found directly from the expansion (I135P 
in the limit (3 — > oo with fixed e<l. Since A — > ir/e 3> 1 in this limit, Eq. (I135P takes the 
form 

w **( e ) 



4vr 3 Af 
c 



Lis (e~" /e ) - (5e/4vr) Li^e - ^) + • • • 
47r3A4 e 5 / 2 ex P {(-7r/e)[l+(5/47r 2 )e 2 + 0(e 3 )]} , e«l. (171) 



This differs from the result (11701) by the small pre-exponential factor yfe and also by small 
subleading terms. The Schwinger formula (I159P is therefore underestimated by the locally 
constant-field expression (I165P for small e as it shown in Fig. 6a. Otherwise, for ne < 1, the 
expression under the sum sign in the first line of Eq. (II 71 ft can be represented as 



£ 5/2 



5 105 , 

e + — = = e 2 



Ann 32n 2 n 2 



e 5/2 



n 5/2 

where 

115 5 



exp (— nn/e) ~ exp [— mr/e*(n)] , ne 1 , (172) 



n" 



e\ ne < 1 , (173) 



e*(n) e An 2 n 2 27r 3 n 3 

can be interpreted as a mass shift via the coherent production of n pairs of particles p2] • 
For strong electric fields e ^> 1, the asymptotic behavior of the rate ( 11591) can be found 



most easily from the integral representation ( 11681) . To this end, we separate the last integral in 
Eq. ( I168P further into three parts 

Li 2 (e^)=/ M-g— - dO^—-- dQ-—. (174) 



& - 1 ./o e e - 1 e - 1 
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Here the first integral is the number T(2) £(2) = 7r 2 /6, the second is represented for 9 < n/e <C 1 
by the expansion n/e — 7T 2 /4e 2 + ■ ■ • ~ n/e, the third is equal to — ln(l — e _7r//e ) and can be 
expanded further as — ln7r/e + 7r/2e + ■ ■ • . Collecting all terms together yields the expansion 

7T 2 

6 eve 

With Eq. (I175j) the large-e expansion of the constant-field rate (11591) reads finally. 

^2 



Li 2 (e 



-7r/e\ 



-6n--l)+G(e- 2 ) , e»l. (175) 



c 



4tt 3 A! 



e 2 Li 2 (e"^) 



c £2 



4vr 3 Af 



7T" 7T /, 7T 

— + - ln--l 

6 eve 



0(e- 2 ) 



e > 1. 



;i76) 



The large-e behavior of the rate fl 1 6 5 [) we find from the expansion (I145P in the limit — > oo 
with fixed e ^> 1. In this limit, substituting A — > n/e <C 1 into Eq. (j!45p yields 



w lcf (e) ~ 



+ 



Ti- 



-In - - 0.5794415 . . . 
2 e 



— fin - + 7 ) - 0.5411802 . 
16 V e 



+ o 



e > 1, 



(177) 



where 7 ~ 0.5772157 ... is the Euler constant |58j. The large-e asymptotic fl 1 77j) is compared to 
the corresponding asymptotic (11761) in Fig. 6b. We observe that the Schwinger formula (11591) 
is overestimated by the locally constant-field expression (11651) as long as e ^> 1, i.e. in the 
strong- field limit. 



0. 00035, 1 1 1 1 -1 1 1 1 -1 0.9 




Figure 6: The dimensionless probability w ■ (Ag/c) to produce an e + e~ pair within one Compton 
space-time volume X\jc m 10 _58 m 3 s as a function of the constant electric field e = Eq/E c . The 
weak- field regime e 1 is shown in (a), the strong-field e ^> 1 in (b). The solid red and blue 
lines refer to the Schwinger formula (I159p and the expression (I165p . respectively. 
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